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1.  INTRODUCTION 

A  workshop  on  the  scientific  results  of  the  ONR-sponsored  Gulf  of  Cadiz  Expedi¬ 
tion  was  held  April  9-11,  1991,  at  the  Applied  Physics  Laboratory,  University  of  Wash¬ 
ington.  The  objectives  of  the  Gulf  of  Cadiz  Expedition  were  to  observe  the  vortices  shed 
in  the  wake  of  Ampere  Seamount,  to  survey  Meddies  formed  by  the  Mediterranean 
outflow  near  Cape  St.  Vincent,  Portugal,  and  to  study  the  structure  and  dynamics  of  the 
outflow  plume  west  of  the  Strait  of  Gibraltar.  The  operational  areas  for  the  expedition  are 
shown  in  Figure  1.1.  During  the  Expedition  (September  4-28,  1988),  velocity,  hydro- 
graphic,  and  dissipation  measurements  were  made  using  a  variety  of  instrument  systems 
(XCP,  XDP,  XBT,  XSV,  CTD,  and  ADCP).  Table  1.1  lists  the  data  types,  number  of 
casts,  drops,  or  observations  made,  the  percentage  of  good  drops  (expendables  only),  and 
the  reports  or  publications  describing  the  processing  of  each  particular  data  type. 

Investigators  from  APL-UW,  WHOI,  JHU/CBI,  the  University  of  Lisbon,  and  the 
Spanish  Institute  of  Oceanography  participated  in  the  expedition.  A  workshop  to  discuss 
preliminary  results  of  the  expedition  was  held  at  the  Woods  Hole  Oceanographic  Institu¬ 
tion  in  October  1989.  The  agenda  for  the  present  workshop  is  given  in  Appendix  A. 
Tabl  ^  1.2  lists  present  publications  according  to  research  objective,  and  Appendix  B  con¬ 
tains  a  complete  list  of  contributions  resulting  from  the  expedition.  Participants  in  the 
April  1991  workshop  are  given  in  Appendix  C.  Appendix  D  lists  past  and  future 
workshops. 

The  following  sections  will  present  highlights  of  the  presentations  given  at  the 
workshop.  Results  and  plans  for  related  work  are  also  included. 

2.  AMPERE  SEAMOUNT  COMPONENT 

2.1  Measurements  of  Ertel  Vorticity  Finestructure  near  Ampere  Seamount 

Recently,  it  has  been  suggested  that  the  oceanic  fine  scale  is  dominated  by  stratified 
2-D  turbulence  (vortical  mode)  rather  than  by  internal  gravity  waves  as  traditionally 
thought;  while  large-scale  flows  can  be  distinguished  on  the  basis  of  frequency,  high- 
wavenumber  structure  is  readily  Doppler-shifted  so  its  dynamics  cannot  be  readily  dis¬ 
cerned  from  simple  time  series.  The  Ertel  vorticity  can  be  used  to  distinguish  between 
these  two  balances  since  anomalies  of  this  quantity  can  be  associated  with  stratified  2-D 
turbulence  but  not  with  internal  waves.  The  Ertel  vorticity  can  be  approximated  as 

ll  =  {f+  VxV)-V(B  +  b)=fN^  +  +  b^v,-byU, 

—  2 

and  is  comprise^  of  a  background  fN  ,  two  linear  anomaly  terms  (stretching  fb^  and  rela- 
tive  vorticity  ^  ),  and  three  nonlinear  anomalies.  For  linear  internal  waves,  the  stretch¬ 
ing  and  relative  vorticity  terms  cancel.  Another  useful  distinguishing  characteristic  is 
that  the  vortical  mode  has  higher  relative  vorticity  than  horizontal  convergence;  near- 
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Table  1.1.  Gulf  of  Cadiz  Data. 


Overview 

Kennelly,  Dunlap.  Sanford,  Kunze,  Prater,  and  Drever,  1989;  The  Gulf  of 
Cadiz  Expedition:  RA'  Oceanus  Cruise  202.  APL-UW  TR  8914. 

Ambar,  Bower,  Cantos-Figuerola,  Kennelly,  Kunze,  Lueck,  Baringer, 
Parrilla,  Prater,  Price,  Richardson,  Sanford,  and  Tokos,  1990; 

Report  on  the  Gulf  of  Cadiz  Expedition  Workshop,  OcL  1  &  2, 1989. 
APL-UW  TM  9-90. 

Drifter 

7  tracks 

Kennelly,  Dunlap,  Sanford,  Kunze,  Prater,  and  Drever,  1989;  The  Gulf 
of  Cadiz  Expedition;  RA^  Oceanus  Cruise  202.  APL-UW  TR  8914. 

Drever  and  Kennelly,  1991: 

Simple  Low-Cost  Drogued  Drifter 

To  be  submitted  to  J.  Atmos.  Oceanic  Technol. 

XCP 

184 

87% 

Kennelly,  Prater,  Dunlap,  Kunze,  and  Sanford,  1989: 

XCP  Data  from  the  Gulf  of  Cadiz  Expedition:  RA' 

APL-UW  TR  8925. 

Prater,  1991;  A  Method  for  Depth  and  Temperature  Correction 
of  Expendable  Probes,  7.  Atmos.  Oceanic  Technol.  (accepted). 

XBT 

T-5 

182 

82% 

Kennelly,  Prater,  and  Sanford,  1989:  XBT  and  XSV  Data  from  the  Gulf 

T-6 

7 

100% 

of  Cadiz  Expedition;  RA  Oceanus  Cruise  202. 

T-7 

40 

90% 

APL-UW  TR  8920. 

xsv 

XSV-02 

52 

79% 

Kennelly,  Prater,  and  Sanford,  1989;  XBT  and  XSV  Data  from  the  Gulf 

XSV-03 

3 

100% 

of  Cadiz  Expedition;  RA  Oceanus  Cruise  202. 

APL-UW  TR  8920. 

XDP 

61 

Lynch  and  Lueck,  1989;  Expendable  Dissipation  Profiler  (XDP) 

Data  from  the  Mediterranean  Out-Flow  Experiment:  RA 

Oceanus  Cruise  202  Leg  V.  JHU-CBl  TR  89-01. 

CTD 

148 

Kennelly,  Sanford,  and  Lehman,  1989;  CTD  Data  from  the 

Gulf  of  Cadiz  Expedition;  RA  Oceanus  Cruise  202. 

APL-UW  TR  8917. 

ADCP 

continuous 

Bathymetry 

continuous 

Digitized  and  merged  with  navigation  database. 

SAIL 

continuous 

SST,  SSC,  wind  speed,  and  direction. 
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Table  1.2.  Gulf  of  Cadiz  Program  Publications  and  Presentations. 


OVERVIEW 

Kennelly,  M.A.,  J.H.  Dunlap,  T.B.  Sanford,  E.L.,  Kunze,  M.D.  Prater,  and  R.G.  Drever,  1989;  The  Gulf  of  Cadiz 
Expedition:  RA'  Oceanus  Cruise  202.  APL-LTW  TR  8914,  Applied  Physics  LaboraicH7,  University  of  Washington, 
Seattle,  WA,  1 15  pp. 

Ambar,  I.,  A.  Bower,  A.  Cantos-Figuerola,  M.A.  Kennelly,  E.L.  Kunze,  R.  Lueck,  M.T.  O’Neil  Baringer,  G.  Parrilla, 
M.D.  Prater,  J.F.  Price,  P.L.  Richardson,  T.B.  Sanford,  and  K.S.  Tokos,  1990:  Report  on  the  Gulf  of  Cadiz  Expedi¬ 
tion  Workshop,  Oct  1  &  2,  1989.  APL-UW  TM  9-90,  Applied  Physics  Laboratory,  University  of  Washington,  Seat- 
Ue,  WA,  58  pp. 


AMPERE  SEAMOUNT  COMPONENT 

Kunze,  E.,  and  T.B.  Sanford,  1990:  Vorticity  and  horizontal  divergence  finestructure  near  a  seamount  (abstract). 
EOS  Trans.  AGU.  72(43),1419. 

Kunze,  E.,  and  T.B.  Sanford,  1991:  Measurements  of  ertel  vorticity  finestructure  in  the  eastern  North  Atlantic.  To 
appear  in  Dynamics  of  Oceanic  Internal  Gravity  Waves,  Proceedings,  'Aha  Huliko'a  Hawaiian  Winter  Workshop, 
Peter  Muller,  Ed.,  Hawaii  Institute  of  Geophysics. 

Kunze,  E.,  M.A.  Keruielly,  and  T.B.  Sanford,  1991:  The  depth  dependence  of  shear  finestructure  off  Point  Arena 
and  over  Pioneer  Seamount.  To  appear  in  J.  Phys.  Oceanogr.,  21. 


MEDDY  COMPONENT 

Prater,  M.D.,  and  T.B.  Sanford,  1989:  Observations  of  an  intense,  young  Meddy  (abstract),  EOS  Trans.  AGU, 
70(43),  1159. 

Prater,  M.D.,  and  T.B.  Sanford,  1990;  Generation  of  Meddies  off  Cape  St.  Vincent,  Portugal  (absuact),  Eos  Trans. 
AGU. 71, 1416. 


MEDITERRANEAN  OUTFLOW  COMPONENT 

Price,  J.F.,  and  M.T.  O’Neil  [Baringer),  1988:  The  descending  Mediterranean  plume  (abstract),  EOS  Trans.  AGU, 
69(44),  1265. 

Ambar,  I.,  and  J.  Verrall,  1990:  Dynamics  of  the  Mediterranean  water  plume  between  the  Suait  of  Gibraltar  and 
Cape  St.  Vincent,  Proceedings  of  the  7th  National  Conference  on  Physics  (Portugal),  p.  457.  72(43),  1419. 

Baringer,  M.O.,  and  J.F.  Price,  1990:  A  simple  model  of  the  descending  Mediterranean  outflow  plume.  The  Physi¬ 
cal  Oceanography  of  Sea  Straits,  L.J.  Pratt,  Ed.,  Kluwer  Academic  Publishers,  Netherlands,  537-544. 

Price,  J.F.,  and  M.O.  Baringer,  1990:  A  model  of  the  Denmark  Strait  overflow  (abstract).  Annales  Geophysicae, 
XV  General  Assembly,  European  Geophysical  Society,  Copenhagen,  April  23—27, 1990,  p.  224. 

Baringer,  M.O.,  J.F.  Price,  and  T.B.  Sanford,  1990:  Mixing  in  the  Mediterranean  outflow  plume  during  the  Gulf  of 
Cadiz  Experiment  (abstract).  Annales  Geophysicae,  XV  General  Assembly,  European  Geophysical  Society, 
Copenhagen,  April  23—21, 1990,  p.  224, 

Parrilla,  G.,  and  A.  Cantos-Figuerola,  1990:  The  Mediterranean  outflow  plume  during  the  Gulf  of  Cadiz  Experi- 
mcni  fSept  1988)  (abstract).  Annales  Geophysicae.  XV  General  Assembly,  European  Geophysical  Society, 
Copenhagen,  April  23-27,  1990,  p.  224. 
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Figure  1.1.  The  operational  areas  for  the  Gulf  of  Cadiz  Expedition  included  Ampere  Seamount,  the  Cape  St.  Vincent 
region  (boxed),  and  the  continental  slope  west  of  the  Strait  of  GibraltJir  (sections  A-I). 


inertial  waves  have  comparable  vorticity  and  convergence,  and  high-frequency  internal 
waves  have  much  weaker  vorticity  than  convergence. 

To  determine  the  segregation  between  these  two  dynamical  regimes  in  the  ocean, 
two  velocity  and  temperature  profile  (XCP)  surveys  were  conducted  just  east  of  Ampere 
Seamount  in  the  eastern  North  Atlantic.  The  surveys  spanned  horizontal  scales  of 
4-7  km  and  vertical  wavelengths  of  40-400  m  (Burger  numbers  of  0.05-30).  Just  as  on 
basin  scales,  the  existing  fine-scale  Enel  vorticity  signal  is  dominated  by  stretching  (Fig¬ 
ure  2.1.1).  The  dynamic  signal,  as  characterized  by  the  relative  vorticity  and  horizontal 
convergence,  is  consistent  with  the  GM  internal  wave  model  and  is  dominated  by  near- 
inertial  waves  on  the  resolved  scales  (Figure  2.1.2).  The  prominence  of  stretching  indi¬ 
cates  that  Ertel  vorticity  finestructure  is  associated  with  nearly  passive  density  finestruc- 
ture  of  very  low  aspect  ratio.  This  further  suggests  that  it  was  not  formed  recently  by 
flow  separation  at  the  seamount  but  might  be  an  artifact  of  subduction  of  surface  mixed 
layers.  Anomalies  of  subducted  water  would  be  long  lived  because  of  the  weak  turbulent 
mixing  in  the  main  pycnocline. 

Section  2.1  was  contributed  by  Eric  L.  Kunze,  University  of  Washington. 
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Figure  2.1.1.  Vertical  wavenumber  spectra  from  survey  1  with  horizontal  fitting  scale 
Ar  =  4  km  for  the  two  linear  and  three  nonlinear  terms  of  the  Ertel  vorticity  anomaly: 
stretching  relative  vorticity  (v^-UyJN^,  {v^-Uy)b^,  tw-isting  b^v^,  and  twisting 
-byU^.  The  spectra  are  presented  in  both  logdog  (left)  and  variance-preserving  (right) 
forms.  The  stretching  peak  at  60-m  wavelength  contains  an  order  of  magnitude  more 
variance  than  the  other  terms  and  thus  cannot  be  explained  by  internal  waves.  At  lower 
wavenumbers  Ck^  >  100  m),  vorticity  and  stretching  become  comparable. 
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Figure  2.1.2.  Comparison  of  the  observed  vertical  wavenumber  spectra  for  horizontal 
convergence  and  relative  vorticity  from  survey  1  (thick  curves)  with  the  GM  model 
integrated  to  horizontal  wavenumbers  of  k/8.6  km”'  (thin  curves).  This  upper  integra¬ 
tion  bound  is  consistent  with  the  fitting  scale  of  the  horizontal  gradients  which  filters  out 
higher  w’avenumbers.  Like  the  GM  model,  the  observed  convergence  exceeds  vorticity 
at  low'  wavenumbers  and  is  comparable  at  high  wavenumbers.  Also,  the  vorticity  peaks 
at  slightly  higher  w'avenumbers  than  convergence  in  both  measurements  and  model.  This 
suggests  that  the  dynamic  (velocity)  signal  is  dominated  by  internal  waves. 
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3.  MEDDY  COMPONENT 


3.1  Dynamical  Aspects  of  the  Cadiz  Meddy  and  Hypotheses  of  Generation 

One  of  the  primary  purposes  (and  accomplishments)  of  the  Gulf  of  Cadiz  Experi¬ 
ment  was  the  detailed  survey  of  a  Meddy  soon  after  its  formation.  The  decay  of  a  Meddy 
has  recently  been  examined  (Armi  et  al.,  1989;  Hebert  et  al.,  1990),  but  the  generation  of 
a  Meddy,  theorized  to  occur  by  one  of  several  different  mechanisms,  had  not  been 
observed. 


Effort  since  the  cruise  has  been  placed  on  a  thorough  description  of  the  Cadiz 
Meddy’s  characteristics  and  dynamical  properties.  The  Cadiz  Meddy  has  two  distinct 
cores,  an  upper  core  at  1050  m  depth  and  a  second,  lower  core  at  1350  m  depth.  The 
upper  core  is  slightly  warmer  and  relatively  fresher  than  the  lower  core,  although  both 
are  similar  in  properties  to  what  is  commonly  referred  to  as  the  "lower"  core  of  the  Medi¬ 
terranean  outflow.  The  upper  core  is  an  isolated  feature,  with  no  connection,  as  seen  in 
temperature,  to  the  Mediterranean  outflow.  This  core  is  in  near  solid-body  rotation  out  to 
6  km  with  a  velocity  maximum  of  0.26  m  s“'  at  9.5  km,  has  a  core  vorticity  near  -0.95/, 
and  is  in  cyclo-geostrophic  balance.  The  upper  core  is  more  energetic  and  axisymmetric 
than  the  lower  core,  which  does  show  a  connection  to  the  outflow.  The  potential  vorti¬ 
city  (2>  a  conservative  tracer  in  the  absence  of  diapycnal  mixing  and  friction,  was  com¬ 
puted  for  the  Meddy  by 


Q 


-  (/k  +  t)‘ 


Vp  ^  1 
P  8 


(/ + 


^  2'’ 

9v 

2  ■ 

f*T 

dz 

4 

- 

and  is  shown  in  Figure  3.1.1.  The  parcels  of  water  in  the  core,  having  near  zero  Q  owing 
to  low  stratification  and  high  relative  vorticity,  are  prevented  from  mixing  outward  by  the 
large  radial  gradient  in  Q. 

We  assumed  that  the  Meddy,  during  its  formation,  entrained  its  two  cores  simultane¬ 
ously,  and  since  then  no  additional  entrainment  or  mixing  of  the  cores  occurred.  Thus 
the  formation  region  is  limited  to  the  locations  where  those  two  water  types  occur  in 
proximity.  Choosing  the  32.15  and  the  32.30  surfaces  as  the  central  isopycnals  of  the 
upper  and  lower  core,  and  mapping  the  potential  temperature  0i  in  the  Gulf  of  Cadiz  on 
those  isopycnals,  we  found  that  the  canyon  south  of  Portugal  was  the  most  likely  forma¬ 
tion  region  (Figure  3.1.2).  In  this  area  the  isobaths  in  the  Gulf  of  Cadiz  converge,  forcing 
the  various  cores  of  the  Mediterranean  outflow  to  constrict  over  abrupt,  canyon  topogra¬ 
phy.  Formation  in  this  region  is  supported  by  extrapolating  the  position  of  the  Meddy 
backward  in  time  using  its  translation  speed.  This  operation  places  the  Meddy  in  the 
canyon  region  30  days  prior  to  our  survey. 
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Figure  3.1.2.  Chart  of  the  Gulf  of  Cadiz,  showing  12  to  12.5°C  waters  on  the  isopycnal 
corresponding  to  the  upper  Meddy  core  and  11.9  to  12.4°C  waters  on  the  isopycnal 
corresponding  to  the  lower  core.  The  two  contours  meet  in  the  canyon  region  south  of 
Portugal  near  36.5°N  and  8.5°W.  The  arrow  extending  northeast  from  the  Meddy  shows 
the  backward-extrapolated  path  of  the  Meddy  for  the  30  days  prior  to  the  survey.  The 
intersection  of  the  temperature  contours,  the  location  of  the  extrapolated  Meddy  path, 
and  the  convergence  of  the  isobaths  all  support  the  hypothesis  that  the  canyon  region 
south  of  Portugal  was  the  formation  region  of  the  Cadiz  Meddy. 
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The  actual  dynamic  mechanism  that  generated  the  Meddy  has  been  limited  to  three 
possibilities:  (a)  frictional  torques  acting  to  reduce  the  potential  vorticity  of  the  Mediter¬ 
ranean  outflow,  creating  a  strongly  anticyclonic  eddy  after  separation  from  topography 
(D’Asaro,  1988),  (b)  the  instability  of  a  buoyancy-driven  boundary  current  (e.g.,  Griffiths 
and  Linden,  1981),  and  (c)  production  of  eddies  by  intermittent  transport  in  an  otherwise 
stable  coastal  current  (Nof,  1991).  The  classical  method  of  eddy  generation  by  geos- 
trophic  adjustment  is  not  considered  to  be  a  viable  mechanism,  since  the  Cadiz  Meddy 
has  an  energy  Burger  number  (Bf  =  KEIAPE)  of  4.5  and  shows  little  evidence  of  having 
been  gravitationally  compressed.  Eddies  generated  by  geostrophic  adjustment  having 
such  a  high  &£  are  compressed  by  a  factor  of  3  or  greater  (McWilliams,  1988). 

Section  3.1  was  contributed  by  Mark  D.  Prater, 

Applied  Physics  Laboratory,  University  of  Washington. 
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4.  MEDITERRANEAN  OUTFLOW  COMPONENT 

4.1  Large-Scale  Structure  of  the  Mediterranean  Outflow 

The  Gulf  of  Cadiz  Expedition  in  September  1988  aboard  the  R/V  Oceanus  partly 
focused  on  the  dynamics  and  large-scale  structure  of  the  Mediterranean  outflow  where  it 
descends  into  the  Gulf  of  Cadiz.  Numerous  closely  spaced  salinity,  temperature,  and 
current  profiles  (acquired  along  sections  A  through  I  in  Figure  1.1)  gave  us  a  greatly 
improved  picture  of  the  overflow,  in  particular  a  detailed  picture  of  the  mixing  and 
dynamics  within  the  first  40  km,  or  500  m  of  descent. 

The  overall  structure  of  the  outflow  is  similar  to  that  described  in  the  historical 
literature  such  as  Heezen  and  Johnson  (1969)  or  Madelain  (1970),  confirming  the  fairly 
steady  nature  of  the  flow  over  decades.  The  outflow  lies  entirely  on  the  continental  slope 
during  its  initial  descent  into  the  Gulf  of  Cadiz,  where  it  mixes  with  fresh  North  Atlantic 
water,  reducing  its  salinity  anomaly  and  increasing  its  transport.  During  the  first  100  km 
of  descent,  the  overflow  transport  doubles  from  0.85  Sv  to  1.9  Sv,  and  the  overflow 
broadens  from  10  km  to  90  km.  Cross-stream  variations  in  the  overflow  T/S  properties 
increase  as  the  flow  spreads.  Slightly  differentiated  water  types  exit  the  strait,  with  the 
saltiest,  coldest  outflowing  water  to  the  south  and  slightly  fresher  and  warmer  outflow  to 
the  north  (Figure  4.1.1).  Since  the  northern,  near-shelf  flow  is  higher  in  the  water 
column,  it  mixes  with  warmer  North  Atlantic  water  than  the  deeper  offshore  flow.  What 
begins  as  less  than  a  0.5°C  cross-stream  variation  in  water  types  in  the  strait  becomes 
more  than  a  2°C  variation  as  the  flow  spreads  within  the  first  100  km  (Figure  4.1.2).  The 
flow  eventually  settles  along  two  preferred  isopycnals:  27.5  and  27.8  (Zenk  1975). 

Our  hydrographic  data  show  a  maximum  core  salinity  of  38.4  psu  at  the  strait  and  a 
decrease  of  0.6  psu  in  40  km.  Combining  hydrographic  data  with  the  current  measure¬ 
ments  suggests  a  much  more  intense  mixing  process.  The  velocity-weighted  salinity  of 
the  outflow  begins  with  a  salinity  of  37.7  psu  near  the  strait  and  decreases  over  1  psu 
within  the  first  40  km.  The  difference  between  the  core  salinity  values  and  the  velocity- 
weighted  salinities  is  due  mostly  to  the  stratified  nature  of  the  outflow.  Core  salinity 
values,  which  are  located  within  10  m  of  the  bottom,  do  not  accurately  reflect  the  average 
propenies  of  the  100  m  thick  outflow.  The  overflow  becomes  neutrally  buoyant  near 
Cape  St.  Vincent,  over  300  km  from  the  strait  where  the  core  salinity  is  36.6  psu.  Since 
the  velocity-weighted  salinity  has  already  attained  a  value  close  to  this  only  40  km  from 
the  strait,  most  of  the  mixing  has  taken  place  within  the  first  40  km. 

The  criterion  for  Kelvin-Helmholtz  instability  is  that  the  gradient  Richardson 
number  is  less  than  1/4.  We  found  low  gradient  Richardson  numbers  in  the  well-mixed 
outflowing  layer  where  the  stratification  is  very  weak  and  also  near  the  interface  between 
outflowing  water  and  the  return  circulation  above.  Ai  the  interface  the  stratification  is 
very  high  but  the  velocity  shear  is  so  large  it  could  overcome  the  stability  of  the  buoy¬ 
ancy  forces.  These  low  gradient  Richardson  numbers  suggest  the  flow  is  unstable  to 
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Section  B 
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Figure  4.1.1.  Potential  temperature  versus  salinity  for  Section  B.  Triangles  represent 
eastward  moving  water  as  determined  from  the  XCPs.  Circles  represent  westward  mov¬ 
ing  flow.  The  symbols  are  scaled  by  the  flow  speed  and  plotted  for  every  4  m  of  data. 
Dashed  lines  show  contours  of  potential  density.  This  figure  shows  rapidly  moving  salty 
Mediterranean  water  moving  to  the  west.  Only  two  stations  across  this  section  indicated 
the  pre.sence  of  Mediterranean  outflow.  The  northern  station  is  very  slightly  warmer  and 
fresher  than  the  southern  station. 
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Section  F 
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Figure  4.1.2.  Potential  temperature  versus  salinity  for  Section  F.  At  section  F,  80  km 
downstream  of  section  B,  there  is  more  than  a  2°C  horizontal  variation  in  outflow  tem¬ 
perature.  The  northern  stations  have  become  warmer  as  the  outflow  mixes  with  warm 
North  Atlantic  water. 
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Kelvin-Helmholtz  instability  and  could  be  mixing  through  this  mechanism.  Bulk 
Richardson  numbers  or,  equivalently,  bulk  Froude  numbers  suggest  the  same  thing  (Fig¬ 
ure  4.1.3).  Supercritical  Froude  numbers,  Fr  >  1,  were  found  near  the  strait  coincident 
with  the  vigorous  mixing  observed  through  the  salinity  and  transport  calculations. 
Farther  from  the  strait,  lower  bulk  Froude  numbers  were  observed  in  regions  in  which  we 
believe  less  entrainment  is  taking  place  (Figure  4.1.4). 

CTD76,  XCP45,  SECTION  C 


Bulk  Richardson  Number: 


Rib  = 


jal-H- 

(AU)2 


Froude  Number: 


F  = 


AU 

(g'H)'-'' 


Rib  =  0.67 
F-1.2 


Figure  4.1.3.  Three-dimensional  velocity  vector  diagram  for  XCP  45/CTD  76.  The 
XCP  velocity  versus  depth  is  plotted  for  every  4  m  of  data.  The  corresponding  salinity 
anomaly  and  gradient  Richardson  number  are  also  plotted.  This  figure  shows  salty  Medi¬ 
terranean  water  moving  toward  the  NW  with  low  gradient  Richardson  numbers.  Low 
gradient  Richardson  numhcrs  are  also  found  in  the  interface  and  the  well-mixed  bottom 
layer.  This  figure  demonstrates  how  bulk  Richardson  numbers  and  Froude  numbers  are 
calculated.  The  average  temperature  and  .salinity  of  the  overflow  were  calculated  using 
the  velocity-weighted  properties,  and  the  change  in  shear  was  calculated  from  vertically 
averaging  the  velocity.  The  bulk  Froude  number  for  this  profile  was  supercritical. 


TM6-91 


15 


Gulf  of  Cadiz 


Z 


Figure  4.1.4.  Station  map  showing  selected  bulk  Froude  number  estimates.  Near  the 
strait,  the  bulk  Froude  numbers  were  higher  than  they  were  at  Section  F,  where  we 
believe  less  mixing  is  taking  place. 


Near  the  strait,  advection,  dissipation,  and  the  Coriolis  force  all  make  important 
contributions  to  the  dynamics  of  the  flow.  Farther  downstream,  the  flow  becomes  a 
damped  geostrophic  current.  After  the  flow  exits  the  strait,  it  takes  a  right-hand  turn  into 
the  Gulf  of  Cadiz.  This  deflection  of  the  current  represents  a  nearly  inertial  turn  in  the 
flow,  where  the  advection  terms  are  the  same  order  of  magnitude  as  the  pressure  gra¬ 
dient.  Examination  of  the  downstream  momentum  flux  divergence  suggests  an  average 
stress  of  2  Pa  retarding  the  flow  within  the  first  100  km  of  descent.  Near  the  strait,  the 
stress  may  be  as  high  as  5  Pa,  while  farther  downstream  the  stress  decreases  to  1  Pa  (Fig¬ 
ure  4.1.5).  Most  of  the  potential  energy  released  from  the  flow  is  due  to  the  descent  of 
the  current  down  the  continental  slope.  This  decrease  in  potential  energy  is  dissipated 
instead  of  increasing  the  kinetic  energy  of  the  mean  flow. 

Section  4.]  was  contributed  by  Molly  O.  Baringer, 

Woods  Hole  Oceanographic  Institution. 
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Energy  per  Unit  Width 


Figure  4.1,5.  Energy  balance  terms  normalized  by  the  width  of  the  outflow.  The  exter¬ 
nal  potential  energy  (dot-dashed  line)  decreases  as  the  outflow  slides  down  the  continen¬ 
tal  slope.  The  changes  in  the  internal  potential  energy  (dashed  line)  mostly  reflect  the 
changes  in  outflow  height.  The  kinetic  energy  (dotted  line)  is  an  order  of  magnitude  less 
than  the  potential  energy.  The  potential  energy  is  not  converted  into  kinetic  energy.  This 
loss  of  total  energy  (solid  line)  suggests  the  presence  of  large  stresses  retarding  the  flow, 
lo  obtain  stress  estimates  from  the  horizontally  integrated  momentum  equation,  we  must 
divide  by  the  width  of  the  outflow.  The  total  stress  is  the  downstream  divergence  of  the 
total  energy  (solid  line).  This  calculation  suggests  a  total  stress  of  5  Pa  retarding  the 
outflow  within  the  first  25  km  of  Section  A.  Farther  downstream  the  stress  decreases  to 
1  Pa. 
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4.2  Modeling  Oceanic  Overflows 

There  are  four  major  overflows  that  supply  dense  waters  to  the  thermocline  or  abys¬ 
sal  North  Atlantic  (Figure  4,2.1).  In  order  of  increasing  density  of  the  final  product, 
these  are  (a)  the  Mediterranean  overflow,  which  injects  warm  salty  water  into  the  North 
Atlantic  thermocline,  (b)  the  overflows  across  the  Iceland-Scotland  Ridge,  which  supply 
warm  and  salty  deep  water  to  the  eastern  North  Atlantic,  (c)  the  Denmark  Strait 
overflow,  which  supplies  bottom  water  to  the  northern  North  Atlantic,  and  (d)  the 
Weddell  Sea  overflow,  which  supplies  cold  fresh  bottom  water  to  the  Weddell  Sea  that 
goes  on  to  become  Antarctic  Bottom  Water.  These  overflows  all  begin  in  marginal  seas 
and  then  spill  over  a  ridge  or  sill  before  they  enter  the  North  Atlantic.  Along  the  way 
they  entrain  oceanic  waters  having  quite  different  properties  and  thus  end  up  with  final 
properties  (T,  S,  and  density)  that  are  quite  different  from  those  found  near  the  sill  or 
strait  connecting  with  the  marginal  sea. 

We  have  developed  a  simple  model  of  overflows  that  we  are  using  to  try  to  under¬ 
stand  why  some  overflows  produce  bottom  water  (the  second  two  mentioned  above), 
while  others  form  interflows  in  the  thermocline.  This  is  especially  intriguing  in  the  case 
of  the  Mediterranean  and  the  Denmark  Strait  overflows  since  the  Mediterranean 
overflow  begins  with  a  very  high  initial  density. 

Our  model  suggests  that  the  principal  external  variables  that  determine  the  intensity 
of  mixing  are  the  bottom  slope  (large  slopes  cause  stronger  currents  and  more  mixing), 
the  eanh’s  rotation  (smaller  rotation  gives  larger  geostrophic  currents),  and  the  initial 
transport  (small  initial  transport  leads  to  thinner  overflow  plumes  which  are  more  suscep¬ 
tible  to  mixing).  A  greatly  simplified  model  including  just  these  variables  can  be  used  to 
find  the  explicit  dependence  of  mixing  on  the  geostrophic  Froude  number  and  the  initial 
Froude  number  (Figure  4.2.2).  Given  an  estimate  of  the  geostrophic  Froude  number,  we 
can  then  use  this  result  to  estimate  what  fraction  of  the  initial  density  difference  will  be 
lost  due  to  mixing,  and  we  have  done  this  for  each  of  the  overflows  shown  schematically 
in  Figure  4.2.1.  In  Figure  4.2.3,  we  plot  the  density  lost  against  a  density  ratio  that  mea¬ 
sures  the  stratification  between  the  initial  overflow  and  the  bottom  water.  If  the  density 
loss  is  greater  than  this  density  ratio,  then  the  overflow  will  form  an  interflow  (upper  left 
half  of  the  diagram);  otherwise  it  will  reach  the  bottom.  The  Mediterranean  overflow  is 
solidly  within  the  interflow  half,  whereas  the  Denmark  Strait  case  is  on  the  other  side, 
though  not  by  a  great  deal.  In  any  event,  we  can  see  the  Mediterranean  case  has  both  a 
moderately  large  geostrophic  Froude  number  and  fairly  strongly  stratified  ambient  water, 
which  both  contribute  to  making  an  interflow. 

Section  4.2  was  contributed  by  James  F.  Price, 

Woods  Hole  Oceanographic  Institution. 
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Oceanic  Overflows 


DS  -  Overflow  of  Arctic  Intermediate  Water  through 
80°W  60°  40°  20°  0°  20°  40°  60°E  Denmark  Strait.  Forms  bottom  water  in  northern 

North  Atlantic  and  contributes  to  the  Deep  Western 
Boundary  Current. 

Before  After 

Overflow  Overflow 

0  0.5°C  1 .0-2.0°C 

S  34.92  psu  34.90  psu 

Oe  28.02  27.92  AOe=0.10 

ISR  -  Overflow  of  Norwegian  Sea  Deep  Water 
through  Faeroe  Bank  Channel  and  across  the 
Iceland-Scotland  Ridge.  Forms  bottom  water  in 
northeastern  North  Atlantic  and  contributes  to  the 
Deep  Western  Boundary  Current. 

Before  After 

Overflow  Overflow 

e  0.0°C  2.5-3.5°C 

S  34.92  psu  35.00  psu 

Oe  28.05  27.90  AOe  =  0.15 

M  -  Overflow  of  Levantine  Intermediate  Water 
through  the  Strait  of  Gibraltar.  Forms  a  warm 
salinity  maximum  within  the  North  Atlantic  thermo- 
cline. 

Before  After 

Overflow  Overflow 

0  13.5.0°C  12.0°C 

S  38.2  psu  36.4  psu 

Oe  28.78  27.69  AOe  =  1.1 

WS  -  Overflow  of  Ice  Shelf  Water  from  RIchner 
Depression  which  forms  Antarctic  Bottom  Water. 
When  combined  with  other  Antarctic  overflows, 
this  water  covers  most  of  the  abyssal  Atlantic  and 
world  ocean. 

Before  After 

Overflow  Overflow 

0  -1.9°C  -0.5°C 

S  34.68  psu  34.65  psu 

Oe  27.94  27.86  A09  =  0.08 

Figure  4.2.1.  Schematic  of  the  major  overflows  that  contribute  to  deep  and  intermedi¬ 
ate  waters  in  the  North  Atlantic. 
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Figure  4.2.2,  Fractional  density  loss  of  the  overflow  as  a  function  of  the  geostrophic 
and  initial  Froude  numbers.  Each  curve  represents  a  different  initial  Froude  number. 
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Figure  4.2.3.  Density  loss  versus  density  ratio  for  the  overflows  of  Figure  4.2.1.  Cases 
that  fall  in  the  upper  left  quadrant  are  expected  to  make  interflows,  whereas  those  in  the 
lower  right  are  likely  to  make  bottom  water,  ("storf."  =  Storfjorden  Fjord,  "den  st"  = 
Denmark  Strait;  "fix'"  =  Faeroe  Bank  Channel) 


TM  6-91 


21 


4.3  Bottom  Stress 


The  stresses  impaned  to  the  Mediterranean  v.’2t;v  plume  from  underneath  are  pro¬ 
duced  by  viscous  buitom  stress  and  form  drag.  The  1  itter  contribution  has  been  reported 
by  other  investigators  (Chriss  and  Caldwell,  1982,  on  a  shelf;  Gust  and  Weatherly,  1985, 
for  the  deep  ocean)  and  most  certainly  is  large  here  because  of  the  bottom  roughness 
(Heezen  and  Johnson,  1969).  The  dissipation  measurements  of  the  XDP  (expyendable  dis¬ 
sipation  profiler)  determine  viscous  stress  whereas  the  profile  observations  of  the  XCP 
are  interpreted  in  terms  of  form  drag.  This  experiment  allows  us  to  compare  viscous 
versus  form  drag  stresses. 

The  viscous  bottom  stress,  t^,  has  been  estimated  using  high-resolution  (0.1m) 
profiles  of  the  rate  of  dissipation  of  kinetic  energy,  e,  collected  with  XDPs.  For  a  steady 
and  horizontally  homogeneous  bottom  boundary  layer,  e  decreases  with  increasing 
height,  z,  above  the  bottom  as 


where  u*  =  (tg/p)^^^  is  the  friction  velocity  and  k  =  0.4  is  von  Karman’s  constant 
(Dewey  and  Crawford,  1988).  Nearly  all  profiles  conformed  to  the  functional  form  of 
(1).  The  spatial  distribution  of  Xg  (Figure  4.3.1)  was  coincident  with  the  axis  of  max¬ 
imum  flow  of  Mediterranean  water.  The  stress  increased  downstream  from  site  3  and 
reached  a  maximum  of  2  Pa  (w*  =  0.044  m  s“' )  between  survey  lines  C  and  D  (see  Figure 
1.1  for  line  locations).  This  maximum  occurred  where  the  plume  veered  north  by  90°. 
The  bottom  stress  decreases  farther  downstream,  but  the  sampling  was  inadequate  for  a 
quantitative  description  of  this  trend. 


A  regression  of  u*  against  the  maximum  speed  of  the  plume  gives  a  drag  coefficient 
of  Cp  =  0.9  X  10“^  (Figure  4.3.2),  while  a  regression  against  the  mean  speed  of  the 
plume  gives  Cq  =  2.6  x  10“^  (Figure  4.3.3).  Drag  coefficients  are  traditionally  based  on 
the  speed  at  1  m  above  the  bottom,  but  such  an  estimate  is  not  available  for  the 
Mediterranean  plume.  However,  it  is  clear  that  the  derived  drag  coefficients  consider¬ 
ably  exceed  the  estimate  of  0.6  x  10“^  by  Dewey  and  Crawford  (1988)  using  a  similar 
technique  for  a  tidal  flow  over  the  continental  shelf. 

Using  the  equation  for  the  conservation  of  mean’ kinetic  energy  and  integrating  it 
over  the  height,  H,  of  the  plume,  it  can  be  shown  that 


H 


<£> 
<U  > 


(2) 


is  a  reasonably  good  estimate  of  the  total  stress,  where  <  >  indicates  an  average  over  the 
height  of  the  plume.  This  estimate  of  the  total  stress  is  typically  twice  the  bottom  stress, 
Xfi,  and  ranges  from  1  to  5  Pa  near  the  axis  of  the  plume  (Table  4.3.1). 

The  use  of  XCPs  to  determine  bottom  stress  is  a  new  endeavor,  and  this  experiment 
provides  an  opportunity  to  examine  the  method  and  compare  the  results  with  the  viscous 


22 


TM  6-91 


0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 

U  [m  s“’] 

Figure  4.3.2.  Regression  of  u*  against  the  maximum  speed  of  the  Mediterranean 
outflow  plume. 
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Figure  4.3.3.  Regression  of  u*  against  the  mean  speed  of  the  Mediterranean  outflow. 
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Table  4.3.1  Mediterranean  outflow  stress  estimates. 
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stress  estimates  obtained  from  the  XDP.  The  procedure  is  to  fit  the  observed  velocity 
profile  to  a  log  function.  The  mean  velocity  profile  is  expected  to  have  the  following 
form: 


u  =  —  ln(z/2o) ,  (3) 

K 

where  u*  and  K  were  defined  earlier  and  lo  is  a  bottom  roughness  scale.  The  u*  deter¬ 
mined  by  (1)  and  (3)  are  not  necessarily  the  same.  They  will  differ  for  a  number  of  rea¬ 
sons,  including  the  influence  of  form  drag,  which  is  exhibited  only  in  the  profile  used  in 
(3). 

One  of  the  first  questions  we  asked  about  using  the  profile  method  on  these  data  was 
the  reproducibility  of  u*  over  a  time  series  of  profiles.  There  were  two  sites  that  were 
repeatedly  sampled.  Site  5,  also  denoted  as  station  4  on  section  C,  was  sampled  a 
number  of  times.  Four  drops  were  made  here  early  in  the  experiment  to  determine  if 
there  were  strong  tidal  fluctuations,  and  two  more  were  taken  later  on  section  C  and 
section  I.  The  observations,  spanning  6.5  days,  are  shown  in  Figure  4.3.4.  The  upper  two 
panels  show  the  filtered,  total  water  column  profile.  The  third  panel  shows  the  down¬ 
stream  velocity  component  at  the  finest  vertical  resolution  possible  (a  value  for  every 
XCP  rotation,  about  0.3  m)  plotted  as  a  function  of  the  log  of  height  above  the  bottom. 
The  flow  estimates  are  fitted  to  a  linear  function  by  a  least-squares  method  to  obtain  the 
stress  estimates,  shown  in  the  bottom  panel  expressed  both  as  stress  (t)  and  friction  velo¬ 
city  (w.)-  These  estimates  are  quite  steady:  T  =  5.2±  1.9  Pa  and  u*  =  7.0±  1.2  cm  s“'. 

It  was  expected  that  the  u*  values  obtained  from  the  XCP  and  XDP  observations 
would  be  different.  Dewey  and  Crawford  (1988)  found  that  u*  values  obtained  from  fits 
to  velocity  log  layers  were  consistently  4.5  cm  s“'  larger  than  those  determined  from  dis¬ 
sipation  profiles.  We  observe  a  similar  result  from  the  fits  to  XCP  profiles  and  the  XDP 
dissipation  profiles.  The  two  instruments  were  often  dropped  several  minutes  and  several 
hundred  meters  apart.  The  comparison,  shown  in  Figure  4.3.5,  has  a  slope  of  3.2. 

As  shown  in  Figure  4.3.6,  there  was  a  large  change  in  bottom  stress  as  the  plume 
accelerated  between  sections  A  and  C.  In  section  A  the  plume  is  thicker  and  slower  than 
in  section  C.  Bottom  stress  x  increases  by  about  an  order  of  magnitude  in  the  maximum 
salinity  part  of  the  plume. 

A  summary  of  the  behavior  of  bottom  stress  and  friction  velocity  as  a  function  of 
downstream  distance  is  shown  in  Figure  4.3.7.  The  sections  have  been  aligned  so  that  the 
position  of  the  salinity  maximum  is  taken  as  the  origin.  Most  of  the  action  occurs  around 
sections  C  and  D  in  the  highly  saline  water.  After  section  D,  the  plume  seems  to  spread 
laterally  and  slowly,  probably  forming  the  two  distinct  cores  that  are  characteristic  of  the 
Mediterranean  outflow. 

Section  4.3  was  contributed  by  Rolf  Lueck.  University  of  Victoria, 

and  Thomas  B.  Sanford.  Applied  Physics  Laboratory,  University  of  Washington. 
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Cadiz  Site  5  /  C4 


Figure  4.3.4. 

Velocity  profiles  at  site  5  (a.k.a. 
section  C  station  4)  expressed 
in  magnetic  north  (v)  and  east 
(u)  velocity  components  (top 
two  panels)  for  the  whole  water 
column  and  the  downstream 
component  expressed  as  a  func¬ 
tion  of  height  above  the  bottom 
(third  panel).  Least-square  fits 
are  shown  for  the  lower  10  m  of 
each  downstream-component 
profile  versus  log  height  off  the 
bottom.  The  lowest  panel 
displays  the  calculated  bottom 
stress  T  and  friction  velocity  u* . 
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Figure  4.3.7.  Summary  of  bottom  stress  observations  from  XCP  profiles  for  sections  A 
through  F.  Note  peaking  of  bottom  stress  around  sections  C  and  D  and  the  broadening  of 
strongly  stressed  flow  along  section  E.  Origins  of  cross-stream  coordinates  are  the  posi¬ 
tions  of  the  maximum  salinity  found  on  each  section. 
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4.4  Comparison  of  Velocities  Derived  from  XCP  and  CTD  Sections 

CTD  and  Expendable  Current  Profiler  (XCP)  data  were  collected  across  the  Medi¬ 
terranean  outflow  in  the  Gulf  of  Cadiz  and  in  coastal  waters  off  Ponugal  during  the 
Oceanus  202  cruise  in  September  1988.  For  each  section,  velocities  directly  measured  by 
XCPs  are  compared  with  geostrophic  velocities  calculated  from  the  CTD  data.  The  XCP 
velocities  are  determined  by  measuring  the  electric  fields  induced  by  the  ocean’s  motion 
through  the  earth’s  magnetic  field.  We  calculated  the  component  of  the  XCP  velocities 
normal  to  each  section  for  comparison  with  the  geostrophic  currents. 

Since  both  velocity-measurement  i>.<  imiques  are  accurate  only  to  within  a  constant 
offset,  all  velocity  profiles  are  referenced  to  a  level  of  no  motion  at  approximately 
300-500  m  depth  which  corresponds  to  the  interface  between  North  Atlantic  central 
water  and  the  Mediterranean  water,  which  flow  in  opposite  directions. 

The  velocities  are  shown  contoured  on  a  plot  of  depth  versus  distance  across  the 
outflow.  We  compared  contour  plots  of  direct  and  calculated  velocities  for  sections  D,  E, 
and  F  (see  Figure  1.1)  and  lines  4,  8,  and  17  (see  Figure  4.4.1).  Qualitatively,  the  agree¬ 
ment  between  XCP-  and  CTD-derived  velocities  is  good.  In  general,  the  patterns  of  dis¬ 
tribution  of  inflow  and  outflow  are  similar,  although  in  two  sections  (section  D  and 
line  8)  the  agreement  was  not  as  good.  In  section  D  (Figure  4.4.2),  the  curvature  of  the 
Mediterranean  outflow  may  introduce  some  deviation  from  a  geostrophic  balance;  in 
line  8  (Figure  4.4.3),  the  presence  of  a  Meddy  at  the  southern  end  of  the  section  (left  on 
the  plot)  is  observed  in  the  XCP  data  but  not  in  the  geostrophic  section,  since  the  flow  is 
in  cyclo-geostrophic  balance. 

Section  4.4  was  contributed  by  Jane  A.  Verrall,  Applied  Physics  Laboratory, 
University  of  Washington,  and  Isabel  Ambar,  University  of  Lisbon. 
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Figure  4.4.1  Survey  pattern  for  Meddy  component. 
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Figure  4.4.2.  EM  velocities  from  XCPs  (top)  and  geostrophic  velocities  from  CTDs 
(bottom)  for  Section  D.  Positive  velocities  indicate  westward  outflow  of  Mediterranean 
water.  The  right  side  of  the  plot  is  closest  to  shore. 
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XCP  Velocities  (cm/s)  Normal  to  Line  0 


Geostrophic  Velocities  (cm/s)  Line  8 


F^ure  4.4.3.  EM  velocities  from  XCPs  (top)  and  geostrophic  velocities  from  (TTDs 
(bottom)  for  line  8.  Positive  velocities  indicate  westward  outflow  of  Mediterranean 
water.  The  nght  side  of  the  plot  is  closest  to  shore. 
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5.  RELATED  WORK 


5.1  Secondary  Circulation  in  the  Faeroe  Bank  Channel  Outflow 

Data  from  three  expendable  current  profilers  (XCPs)  dropped  at  the  center  of  the 
sill  of  the  Faeroe  Bank  Channel  and  a  CTD  station  occupied  there  (Figure  5.1.1)  are  used 
to  examine  the  structure  of  the  northwestward  outflow  of  cold,  fresh,  dense  water  from 
the  Norwegian  Sea  through  the  channel  into  the  Atlantic  Ocean.  The  velocity  profiles 
are  made  barotropic  and  decomposed  into  downstream  (311°T)  and  cross-stream  (41°T) 
components.  Downstream  is  defined  as  the  direction  of  the  difference  between  the  mean 
velocity  from  100-200  m  (in  the  upper  water)  and  that  from  650-750  m  (at  the  nose  of 
the  outflow).  The  resulting,  very  similar,  profiles  are  averaged  as  a  function  of  height 
above  the  bottom  to  obtain  a  mean  profile  (Figure  5.1.2). 

A  mixed  Ekman  and  logarithmic  bottom  boundary  layer  is  present  in  the  bottom  120 
m  (Figure  5.1.2).  A  bottom  stress  of  3.4  ±  0.1  Pa  is  estimated  using  velocity  observations 
from  the  individual  unsmoothed  profiles  in  the  well-mixed  log  layer  (30  m  thick  and 
directed  23  ±  4°  west  of  the  downstream  direction).  This  strong  stress  at  the  bottom 
boundary  creates  an  Ekman  layer  in  the  bottom  120  m,  and  thus  a  secondary  cross-stream 
flow  to  the  southwest  there.  Spin-down  times  for  the  outflow  estimated  from  the  log  and 
Ekman  layer  structures  are  9  and  17  hours  respectively.  The  outflow  takes  about 
26  hours  to  flow  through  the  channel,  or  1.5-3  spin-down  times.  Thus  the  strong  bottom 
friction  must  be  important  in  the  dynamics  of  flow  through  the  channel. 

The  shear  at  the  interface  between  the  outflow  water  and  the  water  above  is 
sufficiently  strong  to  overcome  the  stratification  and  generate  Kelvin-Helmholtz  billows 
at  500-600  m  depth  (where  75%  of  the  values  of  the  gradient  Richardson  number  fall 
between  1/2  and  1/8).  A  cross-stream  flow  of  similar  magnitude  to  that  in  the  bottom 
boundary  layer  but  to  the  northeast  is  found  in  this  high  shear  region  at  the  interface.  If 
persistent,  this  circulation  may  explain  the  pinching  of  the  density  field  observed  at  the 
southwest  channel  wall  in  CTD  sections  across  the  channel.  Cross-stream  flow  toward 
the  southwestern  wall  of  the  channel  in  the  bottom  boundary  layer  would  tend  to  pinch 
isopycnals  there,  while  a  cross-stream  flow  in  the  opposite  direction  at  the  interface  com¬ 
bined  with  the  possible  mixing  there  would  tend  to  spread  the  isopycnals  to  the  northeast. 

Section  5.1  was  contributed  by  Gregory  C.  Johnson, 

Applied  Physics  Laboratory,  University  of  Washington. 
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Figure  5.1.1.  The  bathymetry  of  the  Faeroe  Bank  Channel  (contours  are  in  meters). 
The  location  of  the  CTD  station  and  XCP  drops  at  the  center  of  the  channel  sill  (61.43°N, 
8.34°W)  is  indicated  by  large  dot.  The  arrow  indicates  the  downstream  coordinate  axis 
3 1 1  °T,  which  is  aligned  approximately  with  the  isobaths.  Dark  lines  represent  locations 
of  CTD  sections  across  the  channel. 
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Figure  5.1.2.  The  mean  2-m  smoothed  temperature  (T),  downstream  velocity  (V, 
311°T),  and  cross-stream  velocity  (U,  41°T)  profiles  plotted  against  depth.  The 
stratification  of  the  interface  between  the  warm  upper  layer  and  the  colder  outflow  water 
is  strongest  between  500  and  600  m.  Below  about  750  m  in  the  bottom  boundary  layer 
(BBL),  the  temperature  is  homogeneous.  The  outflow  water  moves  downstream  and  the 
water  above  upstream  with  a  strong  vertical  shear  between  them.  Downstream  velocity 
decreases  below  about  700  m,  in  the  BBL.  Significant  cross-stream  velocities  are 
observed  at  the  depths  of  strong  downstream  shear,  the  interface,  and  the  BBL. 
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5  2  Turbulence  and  Mixing  in  the  Strait  of  Gibraltar 


Two  cruises  in  the  strait  were  undertaken  as  part  of  the  Gibraltar  Experiment,  to 
study  the  influence  of  turbulence  and  mixing  on  the  two-layer  exchange  flow  in  the  strait 
and  to  examine  the  characteristics  of  turbulence  in  the  energetic  strait  regime.  The  sur¬ 
veys  in  the  strait  (with  over  1100  microstructure  profiles)  were  carried  out 
October  12-20,  1985,  and  May  4-16,  1986.  Results  from  the  cruises  (sampling  of  tem¬ 
perature,  conductivity,  and  small-scale  velocity,  as  well  as  ship-mounted  velocity  mea¬ 
surements  and  echosounder  profiling  in  the  upper  275  m)  were  analyzed. 

Mixing  accounts  for  the  changes  in  the  inflowing  layer  along  the  strait.  In  the  west 
end,  the  inflow  layer  is  virtually  pure  Atlantic  water.  In  the  total  inflow  of  1.2  Sv 
observed  at  the  east  end,  0.4  Sv  is  unchanged  Atlantic  water,  and  0.8  Sv  is  half  Atlantic 
and  half  entrained  Mediterranean  water.  The  CTD  and  velocity  data  are  used  to  predict 
turbulent  and  mean  transport  across  isopycnals  due  to  entrainment  between  the  inflow 
and  outflow  layers.  We  use  this  estimate  to  predict  buoyancy  fluxes.  We  predict  an  aver¬ 
age  dissipation  rate  e  of  4  x  10"^  W  kg"'  in  the  middle  layer,  which  agrees  fairly  well 
with  our  observations  but  is  higher  than  the  average  dissipation  rate  observed.  The  high 
dissipation  and  mixing  at  the  Camarinal  Sill  are  nearly  enough  to  account  for  predicted 
buoyancy  fluxes.  The  data  set  shows  the  imponance  of  a  thick  interfacial  layer 
(60-120  m)  and  the  noncoincidence  of  the  velocity  and  density  interfaces.  Because  of 
tidal  flows,  velocity  and  density  structure  are  highly  variable,  complicating  analysis  of 
water  mass  change  within  any  short  sequence  of  profiles. 

The  key  high-dissipation  region  is  the  internal  hydraulic  jump  at  Camarinal  Sill, 
which  was  sampled  repeatedly.  Peak  dissipation  rates  observed  at  the  jump  were 
>10"*Wkg''  in  a  10  m  thick  layer.  (See  Figure  5.2.1  for  an  example  of  a  transit 
through  an  internal  jump.)  Overturns  downstream  of  the  jump  were  often  greater  than 
50  m  thick.  A  second  important  source  of  mixing  is  internal  bores  propagating  eastward, 
which  were  also  sampled  repeatedly.  (Figure  5.2.2  shows  a  burst  which  sampled  an 
internal  bore.  Figure  5.2.3  shows  the  contrast  in  profiles  taken  just  before  and  after  the 
passage  of  the  bore.)  The  dissipation  signal  in  the  bores  decays  fairly  rapidly;  samples 
15  km  and  25  km  east  of  the  sill  are  markedly  different. 

Turbulence  in  the  strait  was  extremely  high,  typically  2  to  4  orders  of  magnitude 
higher  than  in  open-ocean  thermocline  samples.  Dissipation  rates  were  highly  variable 
and  could  not  be  parameterized  in  terms  of  the  shear  observed  in  ADCP  profiles.  Thorpe 
and  Ozmidov  scales  are  nearly  equal. 

Velocity  and  density  data  from  the  east  end  of  the  strait  show  that  the  exchange  was 
different  in  October  than  in  May — it  was  submaximal  in  October  and  maximal  in  May. 

Section  5.2  was  contributed  by  Joel  C.  Wesson. 

Applied  Physics  Laboratory,  University  of  Washington. 
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Figure  5.2.2.  Burst  summary  during  the  passage  of  an  internal  bore  (burst  42).  This 
burst  was  taken  in  Tarifa  Narrows,  about  12  kni  east  of  Camarinal  Sill,  on  May  9,  1986. 
In  the  upper  two  panels  the  density  and  velocity  change  dramatically  when  the  bore 
arrives.  In  the  lower  two  panels  and  e  also  change  substantially. 
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5.3  Kiel  Plans — Lagrangian  Observations  in  the  Mediterranean  Outflow  in  the 

Iberian  Basin 

The  tongue  of  Mediterranean  water  in  the  northeastern  Atlantic  has  been  a  research 
topic  at  the  /nstitut  fiir  Meereskunde  (IfM)  for  two  decades.  In  1988,  we  proposed  to  one 
of  our  funding  agencies,  the  Deutsche  Forschungsgemeinschaft,  to  introduce  a  Lagran¬ 
gian  method  to  study  advection  of,  and  mixing  within,  the  plume  of  Mediterranean  water. 
By  courtesy  of  Tom  Rossby  and  his  team  at  the  University  of  Rhode  Island,  we  had  the 
opportunity  to  transfer  the  RAFOS  technology  to  the  IfM  in  Kiel.  This  Lagrangian 
method  uses  neutrally  buoyant  glass  floats  (1.6  m  long,  1 1  liter  displacement)  to  trace  the 
Mediterranean  water  in  situ.  Acoustic  underwater  communication  enables  RAFOS  floats 
to  collect  travel-time  data  of  coded  signals  from  at  least  two  fixed  sound  sources.  This 
information,  together  with  pressure  and  temperature  data,  provides  the  initial  input  for 
the  calculation  of  three-dimensional  trajectories  of  Mediterranean  water  parcels  (Figure 
5.3.1). 

During  the  past  two  years,  we  have  built  up  a  float  group  in  Kiel.  We  have  learned 
to  build  RAFOS  floats  (Figure  5.3.2),  calibrate,  and  ballast  them,  and  we  have  success¬ 
fully  tested  our  instruments.  Presently  we  are  preparing  for  an  experiment  to  be  carried 
out  from  aboard  our  research  vessel,  Poseidon,  in  the  Iberian  Basin  in  late  May  1991. 
We  plan  to  launch  floats  in  the  undercurrent  off  the  continental  rise  and  in  potential  salt 
lenses  to  be  identified  by  Rolf  Kase  and  his  modeling  group.  They  left  Kiel  for  their  stu¬ 
dies  in  the  Iberian  Basin  in  late  March  1991. 

As  already  suggested  in  the  original  proposal  from  1988,  we  plan  to  continue  our 
work  with  repeated  float  deployments.  We  have  asked  the  Deutsche  Forschungsgemein¬ 
schaft  to  extend  our  funding  through  next  year.  We  plan  repeated  float  seedings  from  a 
location  south  of  the  Portuguese  Algarve  Coast  at  intervals  of  about  10  days.  We  would 
like  to  track  the  upper  and  the  lower  cores  of  Mediterranean  water  from  their  virtual 
source  in  the  outer  Gulf  of  Cadiz,  where  the  tongue  loses  its  contact  with  the  sea  floor. 
This  project  will  be  open  to  constructive  cooperation  if  we  get  the  necessary  funding. 
The  offer  is  especially  attractive  because  the  ensonification  of  the  area  in  the  Iberian 
Basin  could  be  utilized  by  an  unlimited  number  of  floats.  According  to  our  present  plans, 
the  sound  array  may  consist  of  up  to  four  moored  sound  sources,  jointly  supplied  by  IfM 
Kiel  and  Ifremer  Brest. 

Section  5.3  was  contributed  by  Walter  Zenk,  Kathy  L.  Schultz  Tokos,  and  Holger  Konig, 
Institutfur  Meereskunde  an  der  Universitdt  Kiel. 
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Figure  5.3.1.  Example  of  a  RAFOS  float  record  from  the  Iberian  Basin.  On  top  and  in 
the  middle,  we  show  time  series  of  temperature  and  pressure  records  during  a  17-day  test 
occurring  in  fall  1990.  The  trajectory  on  the  bottom  has  been  calculated  from  only  two 
functioning  sound  sources.  L  and  S  denote  the  launch  site  and  the  spot  where  the 
recovered  float  was  first  .seen  by  the  ARGOS  satellite  system.  A  vector-averaged  speed 
on  the  order  of  5  cm  s“'  was  determined  for  the  background  field  of  the  Mediterranean 
water  tongue. 
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Figure  5.3.2.  RAFOS  float  production  in  Kiel.  The  glass  floats  were  assembled,  tested, 
calibrated,  and  ballasted  before  shipment  to  Lisbon.  We  plan  to  launch  them  in  the 
Iberian  Basin  in  May  1991. 
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